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Q\ ; ABSTRACT 

0^ ■ We present interferometric imaging at 33 GHz, with the new superextended configura- 

tion of the Very Small Array (VSA) , of a very deep decrement in the cosmic microwave 
background (CMB) temperature. This decrement is located in the direction of the 

t^J- ■ Corona Borealis supercluster, at a position with no known galaxy clusters, and was 

' discovered by a previous VSA survey (Genova-Santos et a!.). A total area of 3 deg 2 

has now been imaged, with an angular resolution of 7 arcmin and a flux sensitivity of 
5 mjy beam -1 . 

| These observations confirm the presence of this strong and resolved negative spot 

*h . at — 41±5 mjy beam -1 (— 258±29 fiK), with a signal to noise level of 8. This structure 

is also present in the WMAP 5-year data. The temperature of the W-band (94 GHz) 
data at the position of the decrement agrees within 1.2cr n with that observed by the 
VSA at 33 GHz, and within 0.2cr n with the Sunyaev-Zel'dovich (SZ) spectrum. 

Our analyses show that it is a non-Gaussian feature in the CMB at a level of 4.8<r. 
The probability of finding such a deviation or larger in CMB Gaussian simulations is 
only 0.19 per cent. Therefore, an explanation other than primordial CMB is required. 
We have considered the possibility of an SZ effect generated in a diffuse, extended 
warm/hot gas distribution. This hypothesis is especially relevant, as the presence of 
such structures, if confirmed, could provide the location for a significant fraction of 
the missing baryons in the Local Universe. 

Key words: techniques: interferometric - galaxies: clusters: general - cosmic mi- 
crowave background - cosmology: observations. 



1 INTRODUCTION 

The baryon density at z = 0, deriv ed from the total bud get 
over the well-observed components (|Fukugita et al.|[l99Sf ). is 
a factor w2 lower than that at high redshift, inferred through 
independent methods, namely big bang nucleosynthesis 

* E-mail: rgs@mrao.cam.ac.uk 



jBurles et alj l200ll h the Lya forest (|Rauch et all Il997h 
and the cosmic microwave background (CMB) primary 
anisotropies (e.g. Rebolo et al. 2004; Dunkley et al. 2008). 
According to the resu l ts of hydrodynamical simulations 
(ICen fc Ostrikerl 1 19991 : iDave et all 1200 ll : ICen fc Ostriken 
2006) a substantial fraction of these missing baryons could 
be located in the so-called 'warm/hot intergalactic medium' 
(WHIM). This is a very diffuse gas phase, arranged in sheet- 
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like or filamentary structures, with temperatures 10 5 T ^ 
10 7 K, typical baryon overdensities in the range 8pb/(pb) ~ 
10 — 30, and length scales of the order of 10 Mpc. Its inter- 
mediate temperature and low density, in conjunction with 
the presence of many Galactic foregrounds and different ex- 
tragalactic contaminants, make the detection of the WHIM 
rather challenging. Several attempts (some of which have 
produced tentative detections) have been carri ed out, either 
by looking for its possible soft X-ray emission ( Soltan et al.l 
120021 ; iFinoguenov et al.ll2003l; IZappacosta et al 1120051) or by 



identifying UV dNicastro et all |2002| . 120051 ) or soft X-ray 
l|Barcons et al.ll2005h absorption lines in the spectra of more 
distant sources. 

The thermal Suny aev-Zel'dovich (SZ) effect 
l|Sunvaev fc Zeldovich|[l972T ) has been proposed, and in fact 
used, as an alternative tool to search for this hidden matter. 
This is a secondary anisotropy of the CMB due to the 
inverse Compton scattering of its photons by high-energy 
electrons such as those located in the extended atmospheres 
of hot gas (feTc ~ 10 keV) of the richest clusters of 
galaxies. The SZ signal is proportional to the line-of-sight 
integral of the electron density multiplied by the electron 
temperature. This effect is well known and many detections 
have been achieved over the last decade in nearby clusters 
(see e.g. Birkinshaw 1999; Carlstrom, Holder & Reese 2002). 
Structures like superclusters of galaxies, where the WHIM 
is likel y to be located, c ould also build up a significant SZ 
effect (jBirkinshawl 1999). Their lower baryon overdensities 
are compensated for by the l ong path-lengths of the CMB 
hoto ns across them. Indeed. [Hernandez-Monteagudo et alj 
2006) concluded from the results of a numerical sim- 
ulation that ~ 15% of the SZ signal on the sky would 
be generated in this kind of structure. Many studies 
have been carried out to identify statistical detections 
of diffuse intra-supercluster (ISC) gas thro ugh its likely 
SZ e ffect, making use o f the COBE-DMR (iBandav et al 
1996 ) or WMAP dFosalba. Gaztafiaga. fe Castandei 
20031 ; iHernandez-Monteagudo fe Rubino-Martfnl 12004 



Mye rs et al .1120041 ; iHernandez -Monteagudo. Genova-Santos fc 
2004 ) datasets. However, there have been no clear detections 
to date. 

In this context, we selected the Corona Borealis super- 
cluster (CrB-SC) for observations at 33 GHz with the VSA 
extended configuration (Genova-Santos et al. 2005, hereafter 
GS05). These were the first targeted observations search- 
ing for extended SZ in the direction of this kind of struc- 
ture. We found a clear negative feature, which we called 
"decrement H" with a flux density —103 ± 10 mjy beam -1 
(-230±23 fiK) and coordinates RA= 15 h 22 m 11.47 s , Dec.= 
+28°54'06.2" (J2000), in a region with no known clusters. 
Our analyses ruled out, at a level of 99.6%, the explanation 
of this decrement based only on a primary CMB anisotropy. 
This led us to consider an extended SZ effect generated by a 
filamentary structure consisting of warm/hot gas. The pos- 
sibility of an SZ effect from an unknown background cluster 
was also considered, but this is rather unlikely due to the 
large angular size of the decrement, which could only be gen- 
erated by a nearby cluster. The absence of significant X-ray 
emission in this region in the ROSAT XRT/PSPC All-Sky 
Survey map set constraints on the physical characteristics of 
the hypothetical filamentary structure, making it very un- 
likely that it could generate such a big SZ decrement on its 



own. Therefore, the most plausible hypothesis is a combina- 
tion of a negative primary anisotropy and an extended SZ 
effect. 

This work was followed by an observational campaign 
in this re gion with the MITO te lescope at 143, 214, 272 and 
353 GHz (|Battistelli et al.ll2006f ). The results of the analyses 
of the three lowest frequency channels, in conjunction with 
the VSA map at 33 GHz, support the hypothesis of a com- 
bination of a primary CMB anisotropy with an extended 
SZ effect, the relative contribution of the latter component 
being 0.251°'^. The decrement seems to be unresolved by 
the MITO beam (FWHM«16'), and is detected in the three 
lowest frequency channels. 

In the present work we have carried out observations in 
the region of this decrement with the newly superextended 
configuration of the VSA, which has a factor ~2 finer angu- 
lar resolution than the extended configuration used in the 
previous observations. 

In Section 2 we present a description of the VSA in- 
terferometer, emphasizing the differences between the su- 
perextended and the previous configurations. In Section 3 we 
describe the observations and the data reduction procedure 
and in Section 4 we present the final maps. Section 5 includes 
the discussion about the possible origins of the decrement, 
and conclusions are presented in Section 6. 



2 THE SUPEREXTENDED VSA 

The VSA is a 14-element heterodyne interferometer, tun- 
able between 26 and 36 GHz with a bandwidth of 1.5 GHz, 
sited at the Teide Observatory in Tenerife, at an altitude 
of 2400 m. Currently it observes at a central frequency of 
33 GHz, with a system temperature of approximately 25 K. 
A detailed descripti on of the instrument can be found in 
IWatson et ah I (|2003l ). Its main target has been the observa- 
tions of the p rimary anisotr opics of the CMB, aimed at esti- 
Atrim Bananidslaower spectrum (|Scott et~aT]|2003l ; iGrainge et all 



120031 ; iDickinson et al] I l2004l ). although it has also been suc- 
cessfully used to observe the S Z effect in nearby clusters of 
galaxies l|Lancaster et al.ll2005l) . 

The VSA has observed with three different sets of anten- 
nas. The compact array, used from October 2000 to Septem- 
ber 2001, was sensitive in the multipole range £ ~ 150 — 900 
and had antenna apertures providing a primary beam of 
4.°6 FWHM. In October 2001 this was replaced with the ex- 
tended array, sensitive in the multipole range I ~ 300— 1500 
and with a 2.°1 FWHM primary beam. In 2005 the an- 
tenna centres were repositioned and the mirrors replaced to 
provide a l.°2 FWHM primary beam. This superextended 
array covers the multipole range £ ~ 350 — 2450, and pro- 
duces a synthesized beam of ~7 arcmin FWHM. Currently 
it is chiefly dedicated to the observations of primordial fields 
with the objective of investigatin g the power excess discov- 
ered by CBI at£= 2000 - 4000 (|Mason et al.ll2003l ). 

Situated next to the main array are the two 3.7-m 
source-subtractor (SS) antennas for simultaneous monitor- 
ing of radio sources, providing fluxes at 33 GHz on these 
sources, which are then subtracted from the main array data. 
This is a two-element interferometer giving a resolution of 
4 arcmin in a 9 arcmin field. 
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Table 1. Summary of the observations with the VSA. We list the central coordinate of the 6 pointings, the total 
observation time, the integration time and the achieved thermal noise (computed from the map far outside the 
primary beam). 



Pointing 


RA (J2000) 


Dec. (J2000) 


^obs 


Ti nt 


Thermal noise 








(hp) 


(hp) 


(mjy beam ) 


CrB-H 


15 23 00.00 


29 13 30.0 


50 


33 


6.8 


CrB-H-spot 


15 22 11.47 


28 54 06.2 


145 


75 


4.2 


CrB-H-spotN 


15 22 11.47 


29 24 06.2 


25 


8 


18.3 


CrB-H-spotS 


15 22 11.47 


28 24 06.2 


25 


11 


15.0 


CrB-H-spotE 


15 23 56.47 


28 54 06.2 


15 


6 


17.9 


CrB-H-spotW 


15 20 26.47 


28 54 06.2 


15 


6 


23.2 



3 OBSERVATIONS AND DATA REDUCTION 
3.1 Observations 

The observations were carried out mainly during the peri- 
ods: i) October and November 2005; and ii) May and June 
2006. We defined 6 pointings, which are listed in Table [T] 
along with their coordinates. The field CrB-H has the same 
coordinates as the corresponding pointing in the previous 
observations with the extended configuration (GS05). It was 
in this field that the decrement H was detected at the high- 
est level of significance. The pointing CrB-H-spot is centred 
on the position where decrement H was found in those ob- 
servations. With the aim of producing a final larger mosaic, 
we defined four additional pointings with coordinates shifted 
in right ascension and declination. In Table Q] we also quote 
both total observation and integration times (the latter in- 
dicating the amount of data retained after flagging), and the 
thermal noise achieved in each field. Most of the observation 
time was spent on the central fields CrB-H and CrB-H-spot. 
The fractional amount of useful data in the fields CrB-H and 
CrB-H-spot exceeds 50%, but in the four additional fields it 
is considerably lower. This is due to the fact that the latter 
four fields were observed at higher hour angles, where the 
quality of the data is worse, and therefore a larger amount 
of data had to be discarded. The noises are consistent with 
the total integration time in each field. 



3.2 Calibration and data reduction 

The absolute flux calibration of VSA data is determined 
from observations of Jupiter, whose brightness temperature 
is taken from th e WMAP 5-year data: r Jup = 146.6 ±0.7 K 
|Hill et al.l 2008), This flux scale is transferred to other cal- 
ibration sources: Tau A and Cas A. A detail ed description 
of the VSA calibration process is presented in I Watson et al.l 
(2003). The specifications for the superextended configura- 
tion (such as the correction for the fact that Tau A and 
Cas A are partially resolved in the longest baselines) are 
essentially the same as those a dopted for the extende d con- 
figuration, and are explained in lDickinson et all (]2004). As a 
further improvement to these processes for the present con- 
figuration, a secondary phase calibration provided by short 
calibration observations (on Cas A in the observations pre- 
sented in this work), interleaved with the main field obser- 
vations, was introduced. 

The dat a reduction pipeli n e follo wed essentially the one 
described in iDickinson et alj (|2004l ). However, in the su- 



perextended configuration data some further systematics 
have been identified. Some of the daily maps showed stripes. 
This is due to the low quality of the data on the shortest 
baselines at extreme hour angles (around ±2h), detected 
only when stacking the whole dataset in hour angle, but 
not on the individual days. By flagging relevant hour angle 
intervals in these baselines the stripes are removed. 

3.3 Source subtraction 

The source subtraction was carried out similarly to GS05. 
This involves extrapolation to 33 GHz of the fluxes of all the 
sources found in th e NVSS-1.4 GHz (ICondon et al.lll998l ) 
and GB6-4.85 GHz (jGregorv et al.lll996T ) catalogues. In the 
present observations, due to the narrower primary and syn- 
thesized beams, we focus on sources closer than 1.2° from 
the pointing centres and brighter than 18 mjy at 33 GHz. 
We found that all the identified sources had measured fluxes 
above this level, as quoted in Tabled These values were sub- 
tracted from the data, and were obtained, using the same 
estimator as in GS05, from the SS observations performed si- 
multaneously with the extended configuration observational 
campaign (during the periods May 2003-February 2004 and 
August 2004-November 2004). 

This strategy has two main drawbacks. First, as the 
detection limit of the GB6 catalogue is 18 mjy, we are miss- 
ing all the sources with inverted spectra between 4.85 GHz 
and 33 GHz, which may be o f the order of 20% of all the 
sources dWaldram et al]|2003l ). This implies that we could 
be missing around one radio source due to this. Second, we 
are subtracting fluxes that have been derived from SS ob- 
servations carried out two years before the observations of 
the main array. Therefore, variability can be important here, 
and to assess it we have applied a x 2 -test to each source. The 
values of \ 2 quoted in Table [2] have been calculated as 

N ° ba / c o \ 2 

2 \ ^ / Mmeas \ /-.\ 

i— i x * 

where N b s is the total number of observations for each 
source, Si and ASi are the value and the error of each in- 
dividual measurement and Smeaa is the final averaged value. 
This tests the hypothesis that the flux-densi ty time-series 
could be modelled as a constant. Following iBolton et al.l 
(2006), we considered as variable the sources for which the 
probability of obtaining an equal or higher value in a x 2 
distribution with N b s — 1 degrees of freedom is less than 
1 per cent. According to this only the source 1522+2808 is 
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Table 2. List of all the radio sources with measured fluxes at 33 GHz above 18 mjy. We quote the extrapolated 
flux at 33 GHz, and the estimated flux at 33 GHz using the SS observations, which have been used to subtract the 
sources from the data. In the last three columns we list the number of individual observations for each source, the 
fractional variation (per cent) of the fluxes on these individual observations and the result of the \ 2 test. 



Name 


RA 


Dec. 


5"cxt 


Smcas 


Nobs 


AS/S 


x 2 




(J2000) 


(J2000) 


(mjy) 


(mjy) 




(%) 




1516+2918 


15 16 40.9 


29 18 25.0 


37 


37 + 7 


16 


50.4 


6.0 


1519+2746 


15 19 51.0 


27 46 23.0 


33 


34 + 16 


4 


82.1 


3.1 


1521+2830 


15 21 04.8 


28 30 29.0 


87 


34 + 4 


32 


62.5 


36.3 


1522+2808 


15 22 48.9 


28 08 51.0 


135 


97 + 3 


118 


40.0 


172.3 


1523+2836 


15 23 27.9 


28 36 51.0 


62 


32 + 3 


61 


78.9 


45.8 


1524+2900 


15 24 04.2 


29 00 22.0 


57 


32 + 3 


76 


65.2 


45.7 


1527+2855 


15 27 44.7 


28 55 24.0 


38 


35 + 7 


17 


68.8 


12.1 



classified as variable. Its fractional variation, calculated as 

AS 1 / Eff7(ASi) 2 / EMA&f f 
5 mcas sWas V Ei/(A50 2 VEi/(asov ' 

(2) 

is 40%. However, this source lies at a distance of 0.77° from 
the CrB-H-spot pointing (this is the most important point- 
ing in our analysis, as will become clear), where the primary 
beam attenuation is a factor 0.32, so the consequences of its 
variability are ameliorated. To check this we have subtracted 
this radio source from the observed map using the two ex- 
treme fluxes given by the estimated flux plus and minus its 
variability, and found that the flux difference in the centre 
of the map is only 1.3%. 

We have estimated the confusion noise due to unre- 
solved sources below our source subt raction t hresho ld of 
18 mjy by using the formalism of IScheuerl (|l957f ) and 
the best-fitting power-law model obtained from the source 
count of the whole sample of sources in GS05. We found 
o"sour = 3.8 mjy beam -1 - this value is lower than the ther- 
mal noise in any of the fields. 



4 VSA MAPS 

Daily observations, after being calibrated and reduced, are 
held as visibility files. These files are loaded individually into 
MPS, where the final stack for each field is made. Source 
subtraction is then applied in the aperture plane, and maps 
are produced for each pointing using standard AIPS tasks. 
Since we are interested in extended structures, we typically 
use natural weighting. The deconvolution of the synthesized 
beam is carried out using standard AIPS tasks down to a 
depth about the noise level. CLEAN boxes are placed around 
the position of decrement H, and around the radio sources 
in the non-source-subtracted maps. 

The CLEANed maps built up from pointing CrB-H-spot, 
before and after source subtraction, are shown in the top 
panels of Fig. [T] Decrement H is clearly detected (signal-to- 
noise level of 8), with a peak negative flux density —41 ± 

5 mjy beam -1 (with a 7.7 x 6.5 arcmin 2 beam; brightness 
temperature of -258 + 29 fiK) at RA= 15 h 22 m 11.47 s , Dec.= 



+29°00'06.2" (J2000) Q, and with angular size of about 30 x 
20 arcmin 2 . 

In order to enhance the large angular scales, we have 
applied a Gaussian taper function of width a = 150A to 
the visibilities. In the resulting map, which is shown in the 
bottom-left corner of Fig. [T] the detection of decrement H 
is even more remarkable, with a minimum flux density of 
-101 ± 18 mjy beam -1 (with a 15.3 x 10.5 arcmin beam; 
brightness temperature of —197+35 /^K). At this resolution, 
around the peak decrement there are ~ 2 contiguous resolu- 
tion elements each with negative flux exceeding 5cr n (being 
o"n the thermal noise of the map), corresponding to an inte- 
grated decrement flux of —180 ± 19 mjy; another measure 
of the overall decrement is that there are w 4 contiguous 
resolution elements with negative flux exceeding 3<r n , corre- 
sponding to an integrated decrement of —277 ± 26 mjy. 

We have also combined the individual overlapping maps 
of the six independent pointings to produce the final source- 
subtracted mosaic shown in the bottom-right corner of 
Fig. [T] In order to avoid the loss of signal-to-noise ratio, 
and given that the synthesized beams of the six pointings 
have similar shapes and orientations, here we have applied 
the CLEANing to the whole mosaic instead of to the individ- 
ual pointings. To this end, we used the synthesized beam of 
the pointing CrB-H-spot and placed a CLEAN box around 
the decrement. 

In all these maps the decrement is resolved by the 
VSA synthesized beam. This can be clearly seen in Fig. [2j 
where we plot the angular profiles of both the decrement 
and the restoring beam. We have plotted both the radial- 
averaged profiles (left) and the profiles along a line ori- 
ented in northeast-southwest direction (right), where the 
decrement shows the greatest elongation. It can be seen 
that there is some substructure, with a secondary minimum 
towards the northeast, at a position RA= 15 h 22 m 36.64 s , 
Dec - +29°04'06.0" (J2000). 

These results represent a successful confirmation, from 
a completely independent dataset, of the previous detection 
with the VSA extended configuration. In Fig. [3] we show a 

1 This minimum flux value and coordinates (and also other fluxes 
and coordinates presented hereafter) have been obtained from the 
dirty map. This is reliable as the CLEANing procedure only affects 
the regions of the map surrounding the decrement, where the 
CLEAN box has been positioned. Therefore, the minimum flux of 
the decrement may remain the same after the CLEANing. 
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Figure 1. CLEANed maps of the decrement H in the CrB-SC obtained with the VSA superextended configuration. Top: maps of pointing 
CrB-H-spot before (left) and after (right) source subtraction. In the non-source-subtracted map the crosses show the positions of the 
monitored radio sources inside the primary beam FWHM. Bottom left: source-subtracted map of pointing CrB-H-spot, after applying a 
Gaussian taper of a = 150A in the uv plane. The circle indicates the primary beam FWHM (1.2°) of the CrB-Hspot pointing. Bottom 
right: CLEANed VSA mosaic built up from the six pointings listed in Table [T] The noise level is practically uniform across the mosaic at 
a level 5 mjy beam -1 . The units of the colour-scale shown in the top bar are mjy beam -1 . The synthesized beam FWHM is shown 
in the bottom-left corner of each map. 



comparison of the maps in the region of the decrement re- 
sulting from both observational campaigns. Despite that in 
both cases different angular scales are sampled in the aper- 
ture plane of our observations, the two maps are consistent, 
and the shape of the structure is similar in them. In order 
to assess the compatibility of the measurements of the min- 
imum brightness temperature in both observations, we have 
applied a Gaussian taper function of width a = 170A to the 



superextended visibilities. This taper function produces a 
synthesized beam comparable to that of the extended obser- 
vations. The minimum flux decrement in this map is —90 ± 
14 mjy beam -1 (with a 13.9 x 10.3 arcmin 2 beam; brightness 
temperature of —196 ± 30 /iK), whereas the minimum flux 
decrement in the extended map is — 117 ± 15 mjy beam - 
(with a 14.0 x 10.2 arcmin 2 beam; brightness temperature 
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Figure 2. Angular profiles of the decrement H in the CLEANed and source-subtracted maps of pointing CrB-H-spot (solid line) and of 
the restoring beam (dashed line). We represent cases when no tapering is used (top) and when a Gaussian taper of a = 150 A is applied to 
the data (bottom). Left plots show the average flux in rings of increasing radius from the position of the minimum flux of the decrement, 
whereas right plots represent the flux along a NE-SW line. In the left plots the shaded regions show the la error in each ring. 



of -255 ±33 /xK)0. Therefore, both measurements agree to 
within 1.3(J. 



5 ORIGIN OF THE DECREMENT 

In GS05 we presented a detailed discussion about the pos- 
sible origin of decrement H, considering three different sce- 
narios: a primary CMB anisotropy, an unknown cluster of 
galaxies and a concentration of WHIM gas. Our statisti- 
cal analyses yielded a probability of 0.38% for a primary 
anisotropy, while only «0.3 clusters rich enough to produce 
such a deep decrement are expected in the entire surveyed 
region (~ 24 deg 2 ). This cluster hypothesis is further dis- 
favoured by the fact that, due to the absence of any excess of 
X-ray emission in the ROSAT data, a distant cluster would 
be needed, and this is in conflict with the large angular ex- 
tension of the decrement. Furthermore, we have recently car- 
ried out an optical survey with the ING telescope that has 
not shown any evidence of a cluster in the region of the CrB- 
SC (Padilla- Torres et al. 2008, in prep.). A high concentra- 
tion of infrared galaxies at a redshift z ~ 0.11 was found, 
but with a physical distribution unlike that of clusters of 
galaxies. The absence of X-ray emission also set constraints 

2 Note that these fluxes and temperatures have been corrected 
by the primary beam response, which is necessary here as the 
pixel with the minimum flux lies at different distances from the 
pointing centres in the extended and superextended observations. 
This is why the extended flux and temperature values are different 
from those presented in GS05, where no primary beam correction 
was applied and the CLEAN map was used. 



on the physical characteristics that a WHIM structure may 
have. This indeed rules out the possibility of a diffuse struc- 
ture, with a depth of the order of the maximum separation 
along the line of sight between clusters in the core of CrB- 
SC (r;40 Mpc), and with electron temperature typical of 
WHIM. Therefore, we concluded that the most plausible 
hypothesis is a combination of primordial CMB fluctuations 
with an SZ effect from a large-filamentary structure oriented 
in a direction close to the line of sight. We now revisit these 
analyses after bringing in the new superextended configura- 
tion data. 

In order to assess the possible contribution from the 
primordial CMB to this deep decrement, we repeated the 
Monte Carlo analysis described in GS05, but now using as 
template in the aperture plane the visibility points sampled 
by the CrB-H-spot pointing. We carried out 15000 simula- 
tions, adding to each three components: primordial CMB, 
thermal noise and residual radio sources below the sub- 
traction threshold of 18 mjy. For the first of these com- 
ponents, the input p ower spectrum w as generated up to 
t = 3000 with CAMB l|Lewis et al.ll2000l ). from a cosmologi- 
cal model defined by the following parameters: Qb ~ 0.044, 
n M = 0.25, Ov = 0.75, h = 0.73, r = 0.14, W 10 As = 23, 
ns = 0.97, as deri ved from the most recent VSA results 
ilRebolo et ail 120041 ). plus T CM b = 2.725 K (|Mather et ail 
1999) and f2„ = 0. For the thermal noise we performed re- 
alizations using the actual noise values in the visibility tem- 
plate. Finally, the residual radio sources were simulated, in 
the aperture plane and within 1.2° of the field centre, follow- 
ing the source count derived from the flux values obtained 
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Figure 3. Comparison of the source-subtracted and CLEANed 
maps of decrement H obtained with the VSA extended configu- 
ration (contours) and with the VSA superextended configuration 
(colour scale). The units of the superextended map, shown in the 
top bar, are mjy beam -1 . The difference between contours in the 
extended map is 18 mjy beam -1 (1.5cr n ). The circle around the 
centre and the ellipse in the bottom-left corner indicate respec- 
tively the FWHM of the primary and synthesized beams of the 
superextended map. 
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Figure 4. Histogram comparing the flux distribution inside the 
1.2° FWHM of the primary beam in the 15000 simulations (dot- 
ted curve) and in the real data of pointing CrB-H-spot (solid and 
dashed curves, respectively for the dirty and CLEAN maps). We 
also show the la error bars of the simulations (dashed region). 



by the SS in the sources observed in this region (equation (1) 
ofGS05). 

We performed a fluctuat ion analysis similar to 
iRubmo-Martm fc Sunvae^l (|200ot l. In Fig. 0] we show a log- 
arithmic plot of the histograms, i.e. the P{D) functions, of 
pixel values inside the primary beam 1.2° FWHM from the 



15000 simulations (dotted curve) in comparison with the cor- 
responding distribution of the real data (solid curve). For 
negative fluxes below ~ —20 mjy beam -1 the real data 
shows a clear excess with regards to the simulations, caused 
by the presence of decrement H. In the positive flux tail a 
similar excess is found, which stems from the positive spots 
around the decrement, especially from the one located to- 
wards the southeast (see Fig.H). However, this is an artificial 
effect, as we explained in detail in GS05. On one hand, the 
response of an interferometer has zero mean, so the pres- 
ence of a strong negative feature in the map may enhance 
the neighbouring positive spots. On the other hand, the con- 
volution of the synthesized beam may give rise to positive 
features around the strong negative decrement. Note here 
that we have not applied the deconvolution of the synthe- 
sized beam to the simulations, so in Fig. [4] we have repre- 
sented the flux distribution on the dirty maps. When the 
flux distribution on the CLEAN map is plotted (dashed line) 
the excess of the real data in the positive flux tail clearly 
vanishes, while the excess in the negative flux tail remains. 

The percentage of realizations in which the minimum 
CMB flux value is below that found in the real map is 
0.19%. Although a factor of two lower, this value is of the 
order of that obtained when the equivalent study was ap- 
plied to the extended configuration data (GS05). We com- 
puted in the simulations the standard deviations of all pix- 
els within the primary beam FWHM in order to estimate 
the confusion level introduced by each component, finding: 
ccmb = 6.6 mjy beam -1 , a n = 4.8 mjy beam -1 and 
fsour = 2.5 mjy beam -1 . Adding these in quadrature we 
obtain a = 8.6 mjy beam -1 , from which we conclude that 
decrement H is a 4.8er deviation. 

We have also computed the power spectrum of point- 
ing CrB-H-spot, fo l lowing the procedure described in 
iHobson fc Maisingejj (|2002l ). The result is shown in Fig. [5] 
along with that obtained from the l ast published primor - 
dial CMB observations with the VSA (|Dickinson et al.ll2004l ) 
and, for comparison, the theoretical power spectrum that we 
have used in the simulations. In the first bin, at i « 500, the 
power spectrum of the CrB-H-spot pointing shows a 4. Oct 
deviation from the pure primordial CMB behaviour. This 
is caused by the presence of the decrement, as this angular 
scale corresponds to its size, and is similar to the 2.4cr devi- 
ation at I » 550 obtained from the extended configuration 
data (GS05). 

These analyses clearly highlight the non- Gaussian na- 
ture of this decrement, and show that it is unlikely to be 
caused by a primordial CMB fluctuation, confirming the re- 
sults obtained in GS05. Therefore, as the possibility of an 
unknown galaxy cluster seems implausible (see above), the 
hypothesis of a large filament pointing towards us is rein- 
forced. The conclusions about the physical shape, density 
and temperature that this hypothetical structure may have 
are similar to those put forward in GS05. Essentially, a short 
filament located inside the CrB-SC and with an electron 
temperature within the typical WHIM range (0.01 — 1 keV), 
can not account itself for the total SZ signal, without any 
significant X-ray emission. However, we may consider that 
the total decrement could be a combination of a negative 
primary anisotropy and an SZ effect. In fact, this decre- 
ment is detected in the 214 GHz channel (c lose to 218 GHz, 
the null of the thermal SZ effect) of MITO l|Battistelli et all 



8 R. Genova-Santos et al. 



1.4x10" - 

1.2x10*- 

1.0x10* - 

™ 8.0xl0 3 - 

„ 6.0x1 3 : 
+ ■ 
S" 4.0X10 3 - 

2.0X10 3 - 
DC 



/ 1 K^i 



Primordiol CMB fields 
CrB-Hspot superextended ^ 
ACDM ~ " 
ACDM + Residuol RS < 18 mjy 



Figure 5. Power spectrum computed from the pointing CrB-H-spot (filled circles, blue online), in comparison with that obtained from 
the VSA primordial fields observations (Dickinson et al. 2004; open circles, red online). We also plot a ACDM power spectrum (generated 
as explained in section 5), with and without the contribution from residual radio sources below 18 mjy. The strong deviation at i ss 500 
is due to the presence of decrement H. 



2006), indicating that there might be a primordial compo- 
nent. These MITO observations yield that a relative frac- 
tion / = 0.25 of the total decrement is produced by a ther- 
mal SZ effect. However, up to / « 0.7 could be produced 
by this effect within the typical WHIM temperatures and 
without a significant X-ray emission. If, for instance, we as- 
sume that half (/ = 0.5) of the decrement is due to the 
SZ effect from a short filament with electron temperatures 
T c ~ 0.6 — 0.8 keV, then baryon overdensities of the or- 
der Spb/{pb) ~ 1000 — 1200 would be needed. These values 
are in tension with the results of iV-body galaxy formation 
simulations, which predict that most of the WHIM baryons 
(70% — 80%) lie in baryo n overdensities within the range 
5 — 200 (|Dave et al.ll200ll ). However, there is still a signif- 
icant WHIM frac tion at Sps/ips) ~ 1000 (see Fig. 4 of 
iDave et all l|200ll )). and therefore overdensities of this or- 
der are also possible. If we consider a larger filament, then 
these restrictions are relaxed. We could think of a filament 
connecting the CrB-SC with the region at z ~ 0.11 where 
Padilla- Torres et al. 2008 (in prep.) have found a large con- 
centration of infrared galaxies. It would have a length of 
fa 130 Mpc, hold a gas mass M gas w 10 15 M©, and could 
produce an SZ decrement as deep as the structure detected 
in the map (/ = 1.0) without a detectable X-ray emission. 
However, according to the simulations, it is difficult to find 
filaments of this size. This fact, together with the required 
high elongation along the line of sight, make this structure 
rather unlikely. 

It is worthwhile to examine the WMAP maps in the po- 
sition of decrement H. In the top panels of Fig.|6]wc show the 
WMAP 5-year temperature maps corresponding t o the fre- 
quency bands Q, V and W (|Hinshaw et al.l [20081 ). centred 
on the coordinates of CrB-H-spot pointing. The tempera- 
ture towards the centre of these maps seem to be mostly 
negative, at a level of ~ —150 pK. However, it is clear that 
owing to the coarser angular resolution and lower signal-to- 
noise of the WMAP data a detailed comparison with our 
maps is difficult to perform. In order to obtain a more re- 
alistic comparison we have convolved these maps with the 
VSA synthesized beam. To do this we have multiplied the 
WMAP maps by the VSA primary beam, inverse- Fourier- 
transformed them to the aperture plane, sampled the visi- 



bilities in the uv points of the CrB-H-spot stack, and finally 
Fourier-transformed them to the map plane. As a conse- 
quence of this, the resulting maps, which are shown in the 
bottom panels of Fig. [6] contain only the angular scales sam- 
pled by the overlapping of WMAP and VSA window func- 
tions. In these maps the decrement H appears at a high level 
of significance, especially in the Q and W frequency bands. 

The VSA itv-tapered map shown in the bottom-left 
corner of Fig. Q] has a synthesized beam with FWHM~ 
15.3' x 10.5' (in two orthogonal directions). This is compa- 
rable to the 13.2' FWHM of the WMAP W-band beam, and 
indeed both beams subtend the same solid angle. Therefore, 
a comparison between the temperature values derived from 
both maps is reliable. The minimum flux in the VSA map is 
found at position RA= 15 h 22 m 20.62 s , Dec - +28°59'06.2" 
(J2000) and corresponds to a brightness temperature (RJ) 
of —197 ± 35 pK. The thermodynamic temperature in the 
WMAP map on this position is —151 ± 66 pK, which, at 
the frequency of 94 GHz of the W band, is equivalent to a 
brightness temperature of — 120±53 pK. Therefore, the VSA 
and WMAP temperature values agree at a level of 1.2er n . On 
the other hand, taking into account that the ratio of the SZ 
spectrum for the brightness temperature between 33 GHz 
and 94 GHz is 1.510, the WMAP value would correspond to 
a brightness temperature at 33 GHz of — 182±80 pK. There- 
fore, the VSA and WMAP values are compatible within their 
error bars with an SZ effect. 

There are many other mechanisms producing primary 
or secondary anisotropies which are worth discussing here, 
although they usually lead to CMB fluctuations either 
weaker or on smaller angular scales than those produced 
by the thermal SZ effect. 

Kinetic SZ. It is well known that in galaxy clusters the 
kinetic SZ effect is much lower than the thermal com- 
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Figure 6. WMAP 5-year maps llHinshaw et al,ll2008l) in the region of pointing CrB-H-spot. In the top panels the foreground reduced 
maps of bands Q, V and W (from left to right) are depicted. The maps in the bottom panels have been obtained by multiplying and 
convolving the top maps with the VSA primary and synthesized beams, respectively. The units of the colour-scale are fiK in the top 
panels and mjy beam -1 in the bottom panels. 



ponent, typically an order of magnitude smaller. How- 
ever, at lower electron temperatures, like those of the 
WHI M, both components become comparable. In a re cent 
work I Atrio-Barandela. Miicket fc Genova-Santoa (|2008l ) ar- 
gued that the kinetic SZ effect generated in WHIM struc- 
tures could have an important contribution. They estimate 
that 1% of all WHIM filaments of size ~ 5 Mpc, aligned with 
the line-of-sight, produce a combined thermal and kinetic SZ 
effect that could account for the total flux density observed 
in decrement H. In particular, a filament of size L ~ 5 Mpc 
and electron density n c ps fO 3 m~ 3 could produce a com- 
bined thermal and kinetic SZ signal of ~ —200 ^iK. How- 
ever, a structure like this should produce an X-ray signal 

detectable in the ROSAT data. 

Cosmic texture. Also recentlv lCruz et"al . (2007) a rgued that 
the "cold spot" detected in WMAP data (|Vielva et al.ll2004 ) 
could be caused by a texture, a type of cosmic defect due 
to symmetry-breaking phase transitions in the early Uni- 
verse that may produce positive and negative features on 
the CMB (|Turok fc SpergelHl99oT ). This spot has an angular 
size of « f0°, and a minimum thermodynamic temperature 
of ~ —550 fj,K in Q, V and W bands. The predicted num- 
ber of textures is proportional to 8~ 2 (8 being its angular 
size), so many more spots like this are expected in the sky 
at lower angular scales. Therefore, this phenomenon must 
be regarded as an alternative explanation for decrement H. 
Equation (15) of iTurok fc Spergell (|l990T ) allow us to esti- 
mate the number of textures expected in the whole sky of 
a given angular size. Using that formalism we estimate that 



the number of textures of size between 10 and 40 arcmin 
(the angular scale of decrement H) in the whole area ini- 
tially surveyed in the CrB-SC (24 deg 2 ) is 0.3. In the total 
area surveyed by the VSA, combining the CrB-SC survey 
with the primordial fields observations (82 deg 2 ; Dickinson 
et al. 2004), where no features of the size and amplitude of 
decrement H have been found, the total number expected 
is 1.2. Therefore, it would be likely to find a structure of 
this size in the total area observed by the VSA. However, 
it must be noted that different small-scale processes such as 
photon diffusion could smear out textures at sub-degree an- 
gular scales. Current simulations have not enough angular 
resolution to study such small angular scales, so it is not 
clear whether textures of this size could exist. 
Lensing of the CMB. A massive structure like the CrB-SC 
could lens the CMB photons, leading to distortions on the 
mean CMB temperature. However, these distortions are usu- 
ally small. Even in the richest galaxy clusters, CMB lens- 
ing has a characteristic amplitude of a few microkelvins, al- 
thoug h its angular size can extend out to a fraction of a de - 
gree (fseliak fc Zaldarriagd[200til ; lHolder fc Kosowskvll2004l ) . 
In superclusters like CrB, with much lower densities in the 
inter-cluster medium, and whose member clusters are not 
very rich, even lower signals are expected. Therefore, a tem- 
perature fluctuation of ~ —300 fiK like decrement H can not 
be explained in any way by this effect. 

Integrated Sachs-Wolfe and Rees-Sciama effects. When 
CMB photons cross time-varying potentials they are blue 
(in decreasing potentials) or red-shifted (in increasing po- 
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tentials). This mechanism is known as the Integrated Sach- 
Wolfe (ISW) or the Rees-Sciama (RS) effect, when the den- 
sity fluctuation associated with the time-varying potential 
is evolving in the linear or non-linear regime, respectively. 
These effects build up anisotropies in the CMB in the direc- 
tion of extremely massive structures or voids, typically at 
large angular scales (> 5°) and with amplitudes of the order 
of a few microkelvins (see, e.g., Martinez-Gonzalez & Sanz 
1990; Puchades et al. 2006). With such angular scales and 
amplitudes this effect is not likely t o have a significant con - 
tribution to decrement H. However. iDabrowski et ah! (Il999l ) 
studied the physical properties a galaxy cluster may have to 
produce a combined SZ and RS sign al which coul d acco unt 
for the CMB decrement detected by I Jones et al. (1997) m 
the direction of the quasar pair PC 1643+4631 A&B with- 
out X-ray emission, and found a more significant RS signal. 
According to their calculations a galaxy cluster at z = 1 
with T c ~ 1.3 keV and enclosing a total mass of ~ 10 16 Mq, 
would produce SZ and RS effects respectively of the order 
~ 500 /xK and ~ 250 fiK without a significant X-ray emis- 
sion. But, as we have remarked, a distant cluster is not a 
likely explanation for decrement H owing to its large angu- 
lar size. 



6 CONCLUSIONS 

We have carried out observations at 33 GHz with the 
VSA superextended configuration in the direction of a very 
deep decrement previously detected by observations with 
the VSA extended configuration in the CrB-SC (GS05), 
in a region with no known galaxy clusters and without a 
clear excess of X-ray emission. Our results clearly confirm 
the presence of this strong decrement, with flux density of 
—41 ± 5 mjy beam -1 (—258 ± 29 fiK) and coordinates 
RA= 15 h 22 m 11.47 s , Dec - +29°00'06.2" (J2000). It has 
an angular size of about 30 x 20 arcmin 2 , being clearly re- 
solved by the superextended configuration synthesized beam 
(FWHM«7 arcmin). 

We have confirmed that this deep decrement is detected 
in the WMAP 5-year maps, but at a much lower level of sig- 
nificance as a consequence of WMAP's lower angular resolu- 
tion and signal-to-noise. The temperature at the position of 
the decrement in the W-band map agrees at a level of 1.2a n 
with the VSA value, but both temperatures are consistent 
with a SZ spectrum within their error bars. 

A study based on simulations of VSA observations has 
confirmed its non-Gaussian nature. This study has shown 
that it is a 4.8a fluctuation, where sigma here includes 
the confusion noises introduced by the primordial Gaussian 
CMB, thermal noise and residual radio sources. The proba- 
bility of this decrement being entirely produced by a combi- 
nation of these three components is only 0.19%. Therefore an 
alternative explanation needs to be proposed, and we con- 
sidered the most likely to be the SZ effect. A galaxy cluster, 
respectively rich and close enough to produce such a deep 
and large decrement, would probably have already been de- 
tected by existing optical or X-ray surveys. Therefore, this 
SZ effect could be produced in a less dense and hot struc- 
ture such as the WHIM. In order to produce such a deep 
decrement without a significant X-ray emission a filamen- 
tary structure with an elongation along the line-of-sight of 



about 100 Mpc and a baryon overdensity of 500-800 would 
be required. These numbers lie close to the top limit of what 
is expected in these kind of structures, so a filament like this 
would be rather unlikely. However, if we consider that the 
total decrement could be a combination of both a primor- 
dial CMB fluctuation and an extended SZ effect, then we 
could have a structure with a size, density and temperature 
typical of WHIM producing about the half of the minimum 
flux. 

This study reinforces the conclusions addressed in 
GS05. We are confident that our measurements do show 
an excess decrement with regard to the primordial CMB. 
To explain this decrement we require an SZ effect, probably 
generated in a large-scale WHIM structure. If confirmed, 
this would have been the first direct detection of an SZ ef- 
fect from these structures. This would encourage searches 
for similar structures in other superclusters by other CMB 
experiments. If they are found, then a significant fraction 
of the missing baryons in the Local Universe will have been 
identified. 
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